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Universal correlations in spectra of the lattice QCD Dirac operator 
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Recently, Kalkreuter obtained complete Dirac spectra for SU (2) lattice gauge theory both for staggered fermions 
and for Wilson fermions. The lattice size was as large as 12^. We performed a statistical analysis of these data 
and found that the eigenvalue correlations can be described by the Gaussian Symplectic Ensemble for staggered 
fermions and by the Gaussian Orthogonal Ensemble for Wilson fermions. In both cases long range spectral 
fluctuations are strongly suppressed: the variance of a sequence of levels containing n eigenvalues on average is 
given by T,2{n) ~ 2(log n)//37r^ {f3 is equal to 4 and 1, respectively) instead of T,2{n) — n for a random sequence 
of levels. Our findings are in agreement with the anti-unitary symmetry of the lattice Dirac operator for Nc = 2 
with staggered fermions which differs from Wilson fermions (with the continuum anti- unitary symmetry). For 
Nc = 3, we predict that the eigenvalue correlations are given by the Gaussian Unitary Ensemble. 



1. INTRODUCTION 

The study of spectra of the QCD Dirac oper- 
ator is of great importance for the understand- 
ing of the mechanism of chiral symmetry break- 
ing and its restoration above the critical temper- 
ature Tc- Moreover, the fluctuations of the eigen- 
values provide information on the distribution of 
the fermion determinants, and ultimately, on the 
validity of the quenched approximation. In this 
work, we study the Dirac spectrum for SU{2) 
color both with Kogut-Susskind (KS) and Wil- 
son fermions. 

We consider the eigenvalue problems M 
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where D^^ and 

are the usual KS and Wilson 
lattice QCD Dirac operators. The Dirac matrix 
is tri-diagonalized by CuUum's and Willoughby's 
Lanczos procedure The eigenvalues are then 
obtained by a standard QL algorithm. This im- 
proved Lanczos algorithm makes it possible to ob- 
tain all eigenvalues. The accuracy of the eigen- 
values can then be checked by means of sum rules 
for the sum of the squares of the eigenvalues of 
the lattice Dirac operator 



2. SPECTRA OF COMPLEX SYSTEMS 
AND RANDOM MATRIX THEORY 

A basic assumption in the analysis of spectra 
of complex systems is that the variations of the 
average spectral density, p(A), and the spectral 
fluctuations about this average separate. This al- 
lows us to unfold the spectrum {A^} into {A'^}, 
with average spectral density equal to one, ac- 
cording to AJ. — J^'^dXp{X). We study the dis- 
tribution of the number of unfolded eigenvalues, 
{nfc(ri)}, in a stretch of length n, by means of 
its variance, T.2{n), and its first two cumulants, 
7i(ri) and 72 (n). 

These statistics can be obtained analytically 
for the invariant random matrix ensembles of 
hermitean matrices with probability distribution 
given by P{H) ~ ex.p{-TT W H). The matrix el- 
ements can be either real (Gaussian Orthogonal 
Ensemble (GOE)), complex (Gaussian Unitary 
Ensemble (CUE)) or quaternion real (Gaussian 
Symplectic Ensemble (GSE)). They are charac- 
terized by the Dyson index /?, which is equal to 1, 
2 and 4, respectively. The most notable property 
of random matrix ensembles is the stiffness of the 
spectrum, i.e. T,2{n) ~ (2/7r^/3) log(n), instead of 
n for independently distributed eigenvalues. 

In the past decade, spectra of many complex 
systems have been studied in great detail. The 



2 



main conclusion |^ is that, if the correspond- 
ing classical system is chaotic, the spectral cor- 
relations are given by the random matrix the- 
ory (RMT) with the same (anti-unitary) symme- 
tries as the original Hamiltonian. We also note 
that the average spectral density is generically not 
given by RMT. 

3. SPECTRAL CORRELATIONS OF 
THE DIRAC OPERATOR 

For a lattice size of G'^ x 12 and Wilson fermions 
with /3 = 2.12 and k = 0.15, all eigenvalues were 
obtained for 8 independent gauge field configura- 
tions. A histogram of the spectral density (nor- 
malized to 1) in the region [—0.1,0.1] of each of 
these configurations is shown in figure 1. The 
fluctuations of the Dirac spectrum over the en- 
semble and the statistical fluctuations from bin 
to bin seem to coincide. This so called 'spec- 
tral ergodicity', the equality of ensemble averages 
and spectral averages, is a well known property 
of RMT II . We have also calculated the ensem- 
ble average of the number variance of these spec- 
tra unfolded with the ensemble averaged spectral 
density. Within the limited statistics we found 
complete agreement with the GOE prediction. 
For long level sequences much better statistics can 
be obtained from a spectral average rather than 
an ensemble average which we will use in the re- 
mainder of the paper. Below we will show that 
the spectral average of S2 (n) for Wilson fermions 
is given by the GOE as well. 

For SU{2), the anti-unitary symmetry of the 
KS Dirac operator and the Wilson Dirac oper- 
ator is different. We have [T2K, D^^] — and 
['y5CT2K,Qf^] — 0, where K is the complex 
conjugation operator and C the charge conju- 
gation matrix. Because {T2KY = —1, whereas 
{lf>CT2KY = 1 the KS Dirac matrix can be orga- 
nized into real quaternions, whereas the Wilson 
Dirac matrix is real in an appropriate basis. Con- 
sequently, we expect that the spectral correlations 
of the eigenvalues of are given by the GSE, 
and those of are given by the GOE. In Figs. 
I and I we show the number variance, Yi2{ri), and 
the first two cumulants, 71 (rt) and 72(^^)1 versus 
n. The points are obtained from the unfolded lat- 
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Figure 1. A histogram of the spectral density of 
Wilson fermions for 8 independent configurations. 

tice spectra, with the average integrated spectral 
density obtained by fitting a second order poly- 
nomial to a stretch of 500-1000 eigenvalues. An- 
alytical results are represented by the full (GOE) 
and dotted curves (GSE). The KS results are for a 
12^ lattice with 4 dynamical flavors with a mass 
of ma — 0.05. The results for dynamical Wil- 
son fermions were obtained on a 8^ x 12 lattice. 
The values of (3 and k are shown in the label of 
the figure. Clearly, perfect agreement with the 
random matrix prediction is observed. We also 
studied spectral correlations for weaker coupling 
(3 = 2.8, and stronger coupling (3 = 1.8, and 
both quenched and unquenched field configura- 
tions. The number variance was calculated up to 
n = 100 Q The stationarity of the spectral corre- 
lations was investigated by studying correlations 
of stretches of 200 eigenvalues starting from zero. 
In all cases no deviations from the random matrix 
predictions were found. 

4. DISCUSSION AND CONCLUSIONS 

We have seen that the spectral correlations of 
the KS Dirac operator are given by the GSE, con- 
trary to those of the Wilson Dirac operator, with 
the anti-unitary symmetries of the continuum 
theory which are given by the GOE. One might 
wonder what happens to the spectral correlations 
of D^^ in the continuum limit. In a suitable ba- 
sis, we can we write § D^^ = Z?™"*- -l-a^B. The 
operator aP'B is relevant for the spectral fiuctua- 
tions if its norm is larger than the level spacing of 
£)cont. p^^ fixed volume, the level spacing of Aa in 
the bulk of the spectrum is of order 1/N — /V ^ 
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Figure 2. Spectral correlations for Wilson 
fcrmions. 

which is much less than the norm of the perturb- 
ing operator. Only for the low-lying spectrum 
with A ~ 1/V for a ^ 0, a transition to the GOE 
is possible. Apparently, much larger lattices are 
required to study this transition. 

Finally, for SU{3) color, the Dirac operator 
does not have any additional anti-unitary symme- 
tries, and we predict that the spectral correlations 
for both KS fermions and Wilson fermions are 
given by the GUE. In all cases the spectral fluc- 
tuations are strongly suppressed with respect to 
uncorrelated eigenvalues. To some extent, QCD 
is selfquenching. 

The reported work was partially supported by 
the US DOE grant DE-FG-88ER40388. 



Figure 3. Spectral correlations for KS fermions. 
REFERENCES 

1. T. Kalkreuter, Comp. Phys. Comm. 95 
(1996) 1; Phys. Lett. B276 (1992) 485; Phys. 
Rev. D48 (1993) 

2. J. CuUum and R.A. Willoughby, J. Comp. 
Phys. 44 (1981) 329. 

3. O. Bohigas, M. Giannoni, Lecture notes in 
Physics 209, Springer Verlag 1984, p. 1. 

4. M.A. Halasz and J.J.M. Verbaarschot, Phys. 
Rev. Lett. 74 (1995) 3920; J.J.M. Ver- 
baarschot, this proceedings. 

5. A. Pandey, Ann. Phys. 134 (1981) 119. 

6. M.F.L. Golterman and J. Smit, Nucl. Phys. 
B245 (1984) 61; Y. Luo, |hcp-lat/9604025| . 



